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Electro-Mechanical Vibratory System 


Mario Paz 


The problem of predicting the phase angle of two 
self-synchronized rotors starting from rest, is 
presented in this paper. It is shown that with 
insufficient power the rotors may not reach the 
final operating speed of the motors and stay 
locked at one of the lower natural frequencies 
of the vibrating system, thus producing large 
amplitude and failure of the equipment. 


INTRODUCTION 

Mechanical vibrators are widely used for producing the 
motion required in some manufacturing processes such as 
packing cement mixtures, screening, vibration trans- 
portation, vibration pile drivers, etc. The principal 
component of the mechanical vibratory is an unbalanced mass 
rotated by an electric motor. In some of these 
applications, rectilinear vibratory motion is required. 
This motion may be attained using two synchronized motors 
having eccentric masses rotating in opposite directions. 
Synchronization of the rotating masses means coordinating 
the rotation to maintain a definite phase angle between the 
rotation of the two motors. Such synchronization may be 
accomplished by using kinematic coupling such as toothed or 
chain drivers between the two rotors. However, some years 
ago, it became apparent the synchronization could be main- 
tained without resorting to any kind of kinematic coupling. 
It was observed that the angular velocities and phases of 
the motors automatically kept a definite phase rotation 
without the use of coupling. Huygens [1] first observed a 
case of self-synchronization in the seventeenth century. He 
noticed that two pendulum clocks, regardless of their 
initial phase angle, came into synchronization if they were 
mounted on a common resilient support such as a flexible 
beam. He realized that this phenomenon of self- 
synchronization could be explained by the small im- 
perceptible motion in the beam transmitted by the pendulums. 

At the present time, self-synchronization is finding 
wider applications in the design of vibrating machines [2] . 
The use of self-synchronization makes forced synchronization 
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unnecessary. The condition of self-synchronization, as well 
as the phase angle maintained by the rotors, has been 
investigated initially by Blekman [3] for one directional 
motion, by Paz [4,5] for two directional motions and sub- 
sequently by others [6,7,8]. Paz investigated analytically 
and experimentally the problem of predicting the stable 
phase angle and the direction of the resultant force of any 
two rotors with parallel axis of rotation. He showed that 
the phase angle between two rotors, as well as the condition 
necessary for a stable solution could be predicted with the 
application of Hamilton's Principle. 

These developments gave solutions prescribing the 
conditions for s€*lf-synchronization , but did not focus in 
the initial transient portion of the motion until it reaches 
the final steady-state condition. In this transient state, 
the system may have to transverse some of its natural 
frequencies at which resonance may occur, locking the motion 
at a speed lower than the operating speed of the motors, 
thus, resulting in large amplitudes of motion and failure of 
the equipment. The electro-mechanical analysis of vibrating 
equipment during the initial transient motion, is presented 
in this paper. 

CHARACTERISTICS OF INDUCTION MOTOR 

The torque capacity of an induction motor is primarily 
a function of its rotating speed. Fig.l shows a typical 
torque-speed curve for this type of motor. The figure also 
shows an assumed torque-speed curve required to maintain 
rotation as demanded by the load. In this case, the system 
will operate in the transient state until it reaches the 
steady-state condition at which the torque demanded by the 
load equals the torque of the motor, point P in Fig.l. 



Fig. 1 Speed-Torque Curve For Typical Induction Motor 
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At any speed 0 (rad/sec) , the torque differential A T between 
that produced by the motor and the required by the load is 
available to accelerate the motor. Consequently, applying 
Newton's second law for a rotating mass, this differential 
torque is 


d4 

in which J is the mass moment of inertia of the rotating 
parts of the motor and 0 the angular velocity of the rotor. 

MECHANICAL VIBRATING SYSTEM 


The equations of motions of a mechanical vibrating 
system my be obtained by direct application of Newton's 
Law of motion or alternatively by using an energy method 
such as the application of Lagrange's equation. In any 
case, the equations of motion for a linear system may, in 
general, be expressed in matrix notation as 

[MJCYJ + [C]m + [K] CY3 = CF (t) ] (2) 

where [M] , [C] and [K] are respectively the mass, damping 

and stiffness matrices; and [Y3, ($] [V] the displacement, 
velocity and acceleration vectors of the generalized 
coordinates . 


In the solution of Eq. (2), it is convenient to use 
modal superposition method to transform these equations to a 
system of uncoupled equation of the form 

*Z. + 2 t.w.Z. + o) 2 Z. = P. (t) ( i =1 , 2 . . ,N) (3) 

l fx 1 1 1 1 l 

where Z^ = modal displacement 

9 

= modal velocity 
Z^ = modal acceleration 
ui = natural frequency 
^ = modal damping ratio 

and (t) = modal force 

which is given by 


N 


*!<« = I ♦ji 

j-1 

where 4> ^ is the j-component of i-th eigenvector. 


(4) 


The solution of the uncoupled equation, Eq. (3) during 
the transient motion may be obtained by numerical in- 
tegration in terms of the modal coordinates Z^. 
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The motion in function of the generalized coordinates is 
then given by the transformation 

{Y} = [f][Z3 (5) 

where [<£] is the modal matrix of the system. 



Fig. 2 Mathematical model for the example of a vibrating 
system excited by two self- synchronizing rotors. 


TORQUE DEMANDED BY THE LOAD 

The torque demanded by the load may be determined from 
the power consumed by the exciting forces of the rotors. 
The forces F.,(t) and F 2 (t), developed by the two rotors 
exciting the system are giver* by 

Fjl (t) = m'e 0^ (6) 

F 2 (t) = m'e 0 2 

where m' and e are the eccentric mass and its corresponding 
eccentricity; ©^andf^the angular velocities of the rotors. 

The power demanded by the load may be calculated as the 
product of the rotor forces times the corresponding linear 
velocities of the points of application of these forces. It 
is convenient to express this power separately for the two 
force components in the plane of the rotating forces as 


F 11 (t) = F 1 (t) sin 0 1 (t) 


F 12 (t) = F 1 (t) COS 0 1 (t) 


(7) 
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and 


21 (t) 

- p 2 (ti 

sin 0 

2 (t) 

22 (t) 

II 

ct 

cos 0 

2 ^ 


where F^. is the j -component of the i-th motor; q and 0 the 
angular displacements of the rotors. 1 2 


The power demanded by the load is given by 
p n (t) = F 11 (t) .X 11 (t) 

P 12 (t) = F 12 (t) • ^12 (t) 


(9) 


and 


P 21 (t) = F 21 (t) **21 (t) 

P 22 Ct) = F 22 (t) .X 22 (t) 


( 10 ) 


where Pj. = Power of j force component of the i-th motor 
3 . • = Linear velocity of the force component j of the i-th 
motor . 


Finally, the torque Tj. demanded by the load is obtained by 
simply dividing the power components in Eqs. (9) and (10) by 
the corresponding angular velocity of the rotors. Con- 
sequently, the total torque demanded by the loads is given 
by 

2 

T i. - lipid' 4 ' + p k2 (t > 1 (11 » 

k=l K 

ANGULAR VELOCITY OF THE ROTORS 


The angular velocity for the rotors is determined by 
numerical integration of Eq. (1) which expressed in finite 
differences form may be written as 


A 0<t) 


AT( t ) 
J 


At 


( 12 ) 


where A 0(t) is the incremental angular velocity of the 
rotors attained in time step A t and A T(t) is the torque 
available at accelerate the rotors. As illustrated in 
Fig.l, this torque is given by the difference between the 
torque of the motor T (t) and the torque T ft)demanded by the 
load; thus for each motor, the differential torque is given 
by 

A T (t) = T M (t) - T L (t) (13) 
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The angular velocity for each rotor # is then calculated, 
starting with initial velocity zero (d(0)-0), through the 
recurrence formula 


0 (t+A t) =0 (t) + A0 (t) 


(14) 


The angular displacement of the rotors, 0^, and © 2 is 
calculated as 


G 1 (t+ t) = e^tJ + S^t) x A (t) 
© 2 (t+ t) = © 2 (t)+0 2 (t) x A (t) 


(15) 


starting with the initial displacements of the rotors 
0 1 (t=O) and 0 2 (t=O). 

The phase angle <^-(t) between the rotors is then given 
by the difference between the angular displacements of the 
rotors. Hence: 


*C(t) = 0 1 (t) - & 2 (t) 

NUMERICAL EXAMPLE 

Figure 2 shows the mathematical model of a vibrating 
conveyor, excited by two motors with eccentric masses 
rotating in opposite directions. The differential equations 
of motion for the six degrees of freedom in this model were 
obtained using Lagrange's equation. The numerical values 
for the components of this system are indicated in Table 1. 
Discrete values of the torque-speed function for the motors, 
provided by the manufacture, are given in Table 2. 

A time step At =0.01 sec. was selected and damping in the 
system was assumed equal to 5% of the critical damping in 
each mode. 


TABLE 1 NUMERICAL VALUES VIBRATING SYSTEM OF FIG. 2 


BASE: 


Mass: M^= 31 . 863 ( lb-sec^/in) 

Mass moment of inertia: 

I 1 = 42 , 772 ( lb- sec 2 - in) 

EXCITER: 

Mass: M 2 = 15 . 236 ( lb-sec^/in) 

Mass moment of inertia: 

I 2 = 4,578 (lb-sec -in) 
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Springs costant: 


K, (axial ) = 3k (k = 9 , 400 ( lb/in) 

K 2 (axial )= 9k 

(transverse) = 30k 
K 2 (transverse) = 90k 

Distance : 

a = 9.0 in 
b. = 29.0 in 
hi = 13.25 in 
c 2 = 27.0 in 
cU = 8.0 in 

Eccentric rotor: 

m'e = 0 . 1425 ( lb. sec^ ) 


TABLE 2 SPEED -TORQUE VALUES 


Speed 

(rad/sec) 

■ ' 0 ' 

47.12 

162.3 

184.3 

188.5 

Torque 

(lb-in) 

428.4 

321.3 

616.75 

142.8 

0 



Fig. 3 Motor speed and phase angle of the rotors - Motors 
having full torque. 


The analysis of the system modeled as shown in Fig. 2 
provided the following three non-zero values for the natural 
frequencies : 
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=114.7 rad/sec, o) 2 =195.7 rad/sec and w 3 =783.6 rad/sec. 

The solution output includes the motion of the system and of 
the rotors, as a function of time, during the initial 15 
seconds of the transient state. The solution also includes, 
the torques of the motors and the torques demanded by the 
load as well as the phase angle between the rotors in the 
transient state. 

For the motors having the speed-torque characteristic 
depicted in Fig.l (Numerical values in Table 2), Fig. 3 
shows the plots of the phase angle between rotors and the 
angular speed of the motors during the initial 40 seconds 
motion, while Fig. 4 shews these plots assuming that the 
motors can provide only one-half of the torque indicated in 
Table 2. The initial phase angle was 60 degrees due to the 
physical arrangement of the rotors under the effect of 
gravity. 



Fig. 4 Motor speed and phase angle of the rotors - Motors 
limited to 50% of the full torque. 


DISCUSSION 

It may be observed in Fig. 3 that for the motors having 
the speed-torque characteristics of Fig.l, the motors 
rapidly, in about five seconds reach their operating speed 
in the neighborhood of $ =185 rad/sec and that rotors attain 
self-synchronization with a phase angle <<=0 degrees. On the 
other hand, Fig. 4 shows that the motors do not have 
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sufficient power to reach the operating speed and become 
locked in the neighborhood of resonant speed w = 114.7 
(rad/sec) attaining out of phase self-synchronization of 
the rotors (<X-i8G ff ) . In the first case, the self- 
synchronized machine will perform as intended in the design, 
while in the second case the system will vibrate erratically 
at the resonant condition with large amplitudes resulting in 
failure of the equipment. 
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